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In this survey, we review some of the low energy quantum predictions of General Relativity 
which are independent of details of the yet unknown high-energy completion of the gravitational 
interaction. Such predictions can be extracted using the techniques of effective field theory. 


1. INTRODUCTION 

The study of quantum gravity is a major research field. Over the past two decades, there has been a transition in the 
understanding of this subject that has not yet been fully absorbed by scientists outside of the field. It used to be stated 
that general relativity and quantum mechanics were incompatible. There were many reasons given for this conflict, 
some of which look foolish from the modern perspective. However, a modern view is that general relativity forms 
a quantum effective field theory at low energies. As described below, effective field theory is a standard technique 
to describe quantum effects at low energy where one knows the active degrees of freedom and their interactions. 
The effective field theory allows predictions which are valid at those energies. This does not eliminate the need 
to understand gravity at very high energies where many interesting effects occur. However it is still remarkable 
progress, as we now understand that gravity and quantum mechanics can be compatible at the energies that have 
been experimentally probed. 

The effective field theory treatment allows the separation of quantum effects which follow from known low energy 
physics from those that depend on the ultimate high energy completion of the theory of gravity. The key ingredient 
follows from the uncertainty principle in that high energy effects are very local while those from low energy are non¬ 
local. Indeed there are some results that can be described as “low energy theorems” of quantum gravity. This means 
that they are the outcome of any consistent theory of quantum gravity. The only assumptions of that full quantum 
gravity theory is that it limits to general relativity at low energy. Given how much the combination of gravity and 
quantum mechanics has been maligned in the past, it is remarkable that such universal results for quantum gravity 
can now be formulated. 

In this review, we provide a survey of some of the low energy theorems which have been calculated thus far. Most 
of these are scattering amplitudes, as these are the structures that quantum field theory calculates most readily. In 
the future we hope that equivalent results can be developed for other gravitations settings. 


2. EFFECTIVE FIELD THEORY 

Physics is an experimental science. We only reliably know the degrees of freedom and their interactions which have 
been probed by present experiment. This means that there almost certainly exists new physics to be found beyond the 
present energies. For gravity, this caveat is especially relevant since the Planck scale of Mp = = 1.2 x 10^® GeV, 

where G is the Cavendish constant, seems to be the obvious location for new physics. Experimentally, however, we 
are so far from the Planck scale that we have little hope of uncovering the nature of this new physics in the foreseeable 
future. 

Quantum mechanics does seem to care about this unknown new physics, since in perturbative calculations we are 
instructed to sum over a complete set of intermediate states—at all energies—when making quantum calculations. 
This difficulty is solved, however, by a simple mechanism—the Heisenberg uncertainty principle. The effects of very 
high energy appear as short-distance phenomena to us and thus appear as local terms in a Lagrangian. The coefficients 
of these local terms are the residual manifestation of the high energy physics. Most of these Lagrangians are basically 
irrelevant because they are suppressed by powers of the heavy scale. Thus the unknown physics is reduced to a very 
few parameters which we can either measure (if we do not know the high energy theory) or predict (if we think that 
we do know this theory). 

On the other hand, low energy physics is of a very different character, since light particles can propagate long 
distances and their low energy effects are not local. The distinction between local and non-local is the key to separating 
the physics of low energy from high energy. To be sure, even light particles, when treated in loop diagrams, can have 
effects from very high energies since loops integrate over all energies. These high energy effects are not reliable, as 
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we have not yet experimentally probed that part of the theory. But this is not a problem. Again these high energy 
effects are local, and merely are absorbed into the renormalization of the coefficients of the local effective Lagrangian. 

Effective field theory (EFT) is the technique that takes advantage of these physical properties in order to make 
predictions at low energies which are reliable, because they utilize only the low energy particles and their low energy 
interactions. Useful references on EFT and its applications can be found in Hi- 

As a more technical explanation, consider a general quantum field theory with light and heavy degrees of freedom. 
In quantum field theory it is straightforward to represent what is going on by use of path integral methods. If (j), (<i>) 
represents a light (heavy) field respectively, then the functional integral which characterizes the full quantized theory 
is given by (note we are using h = 1 


W = 


J [d(j)] [d$] exp i 


d)). 


( 1 ) 


Now suppose that we integrate out the heavy degrees of freedom. What is left is a functional integral in terms of a 
non-local “effective” interaction which characterizes the theory in terms of only the light degree of freedom (j) 


W = N [d(j)] exp i / d'^xCef / {4>) 


( 2 ) 


but which includes the virtual effects of the heavy degrees of freedom d) to all orders. 

In this setting, it is possible to produce low energy theorems. This phrase refers to calculations of amplitudes or 
relations between amplitudes that remain valid independent of any modification of the high energy component of 
the theory. This can only happen if such relations depend uniquely on the low energy part of the theory. Any high 
energy theory that is capable of generating these particles and these interactions at low energy must yield identical 
results. Effective field theory allows us to calculate these relations. In the context of this review, the conditions for the 
existence of these low energy theorems is that the ultimate quantum gravity theory must reduce to general relativity 
in four dimensions at low energy. Conventional quantum field theory, expressed through path integral quantization, 
is assumed to apply. There are not really any other ingredients. This combination of general relativity and quantum 
field theory automatically behaves as an effective field theory at the lowest energies, and we can extract low energy 
theorems by the use of effective field theory, as shown below. 


3. LOW ENERGY THEOREMS OF QCD 

In order to understand the use of effective theory in a non-renormalizable theory such as general relativity, it 
is useful to first examine low energy quantum chromodynamics (QCD). In QCD at the lowest energies there exist 
only light pions which are cNnamically active and the interactions of these pions are constrained by the original 
chiral symmetery of QCD[I1-Q. The resulting effective field theory—chiral perturbation theory—has many aspects 
in common with general relativity. Chiral perturbation theory has been exceptionally well studied both theoretically 
and experimentally. We provide a somewhat detailed review here in order to set the stage for a parallel treatment of 
general relativity. 

In order to understand low energy QCD, we begin by introducing the property of “chirality”, defined by the 
operators 


TL.fl = 2(1^75) = 2 ^ ( 3 ) 

which project out left- and right-handed components of the Dirac wavefunction of quarks via 

ipL = riV i’R = with Ip = (4) 

In terms of these chirality states the light {u, d) quark component of the QCD Lagrangian can be written as 

^cfcD =Q{iP - 'm)<l = JP<1L + qRi Pqr - QLmqR - qRiriqL (5) 

where , = ( “ ) is a two component sp.pot and ™ is the a.d ppa* „ass „at„x. We note that 

the limit of vanishing mass 


nud 

^QCD 


m —)-0 


qLi pqL + qRi pqR 


( 6 ) 
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is invariant under independent global left- and right-handed rotations 

3 3 

ql -t exp(i^ Tjaj)qL = Lql, qr exp(i ^ Tj(3j)qR = Rqu (7) 

i=i j=i 


This is SU (2)l {E> SU{2)r invariance or chiral SU (2) SU (2). Continuing to neglect the light quark masses, we see 
that in a chiral symmetric world one might expect six—three left-handed and three right-handed—conserved Noether 
currents 


qLJti^nqL, qRlt^^nqn ( 8 ) 

Equivalently, by taking the sum and difference of the left- and right-handed currents we should have three conserved 
polar vector and three conserved axial vector currents 


A* 


1 


Al 


1 


(9) 


The polar vector symmetry is seen in the particle spectrum and is just isospin symmetry. However, the axial sym¬ 
metry is not observed in a Wigner-Weyl fashion—there are no parity doublets—but rather is spontaneously broken. 
The implication of this breaking is that, via Goldstone’s theorem^, there must exist nearly massless pseudoscalar 
particles— pions—which are approximate Goldstone bosons for the chiral symmetry of QCD. Moreover, the pion 
interactions must obey numerous symmetry restrictions. The simplest way to keep track of these symmetry require¬ 
ments is to write a general effective Lagrangian which obeys this symmetry. This can be accomplished by use of a 
nonlinear function of the pion field, 


U = exp 



transforming as U ^ LUR^ 


( 10 ) 


with L, R being the 2x2 SU{2) matrices in SU{2)r^r respectively which were introduced above. Lagrangians 
constructed using U retain the chiral symmetry and obey the symmetry requirements. The general Lagrangian is 
then constructed in an expansion in the number of derivatives, so that the terms with the fewest derivatives will be 
most important at low energy. 

We infer then that the lowest order SU(2) effective chiral Lagrangian can be written as 

£(2) ^ + f/l'). (11) 

where the superscript 2 indicates that we are working at two-derivative order or one power of chiral symmetry 
breaking— i.e. rn^ oc uiu + rrid- This Lagrangian is also unique —if we expand to order cfA 

= Tr-^x • x -^t ■ ■ d^cf ), (12) 

TT TT TT 

we reproduce the free pion Lagrangian, as required— 

= \9f,(f>-d^^cf>-^mlcf>-cj). (13) 

At higher orders interactions are generated. For example, in the case of pion scattering, expanding order 

(jA we hnd 

4 " + ( 14 ) 

which yields for the on shell tttt T-matrix 

'^ab;cd{_Qa-) Qbi Q_ci Qd) — ^tt) ^ab^bdij^ ^tt) ^ad^bci^^ ^tt)] 

-^TT 

Defining more generally 


'^ab]cd (s, t, u) = A(s, t, u)SabScd + A{t, S, u)6acSbd + A{u, t, s)6adSbc , 


(16) 
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Experimental Lowest Order First Two Orders 


an 

0.220 ± 0.005 

0.16 

0.20 


0.250 ± 0.030 

0.18 

0.26 

an 

- 0.044 ± 0.001 

- 0.045 

- 0.041 

b'i 

- 0.082 ± 0.008 

- 0.089 

- 0.070 

a\ 

0.038 ± 0.002 

0.030 

0.036 

b\ 

0 

0.043 


a2 

(17 ± 3 ) X 10 "^ 

0 

20 X 10 “^ 

di 

( 1.3 ± 3 ) X 10 "^ 

0 

3.5 X 10 -'‘ 


TABLE I: The pion scattering lengths and slopes compared with predictions of chiral symmetry. The last column has been 
taken from ^ 


we can write the chiral prediction in terms of the more conventional isospin language by taking appropriate linear 
combinations 


T°(s,t,t() = 3A(s, t, u) + s, u) + A(u, t, s), 

T^{s,t,u) = A{t, s,u) — A{u,t, s), 

T'^{s,t,u) = A{t,s,u) + A{u,t,s). (17) 

Partial wave amplitudes, projected out via 

1 

= 71“ / d{cos9)Pi{cos9)T^{s,t,u), (18) 

647r 

can be used to identify the associated scattering phase shifts via 

1 . 

Then from the lowest order chiral form Eq. (HU 


A(s, t, u) = 


F3 


we generate predictions for the pion scattering lengths and effective ranges @ 


an = 


.0 ^ 

° 47rF2 


„2 

327rF2 ’ ° 

2 


mt 


ml 


lenFl ’ 


247rE2 


to: 


h‘‘ — _ 

~ SttFS ’ 


( 20 ) 


( 21 ) 


comparison of which with experimental numbers is shown in Table 1. These are the start of some low energy theorems 

of QCD0. 

However, we can do better by considering loop effects. These will generate corrections to the tree amplitudes and will 
bring in imaginary parts for the amplitude that are necessary to satisfy unitarity. We will see that these corrections are 
expressed in an expansion in the energy, such that the tree results are the lowest energy results and one loop results are 
corrections to the tree level. Inclusion of loop effects comes with a price—numerous divergences are introduced and this 
difficulty prevented progress in this field for nearly a decadeQ until a paper by Weinberg suggested the solution^ — 
dealing with such divergences, just as in QED, by introducing phenomenologically determined counterterms into the 
Lagrangian in order to absorb the infinities. We show in the next section how this can be accomplished. 


3.1. Effective Chiral Lagrangian for QCD 

We now apply Weinberg’s suggestion to the effective chiral Lagrangian, Eq. m- As noted above, when loop 
corrections are made to lowest order amplitudes in order to enforce unitarity, divergences inevitably arise. However, 
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there is an important difference from the familiar case of QED in that the form of the divergences is different from their 
lower order counterparts—the theory is nonrenormalizable! The reason for this can be seen from a simple example— 
TTTT scattering. In lowest order there exists a tree level contribution from which is 0{p^/F^) or Oimf^lFff) where 

p represents some generic external energy-momentum. The fact that p appears to the second power is due to the 
feature that its origin is the two-derivative Lagrangian £ 2 - Now suppose that tttt scattering is examined at one loop 
order. Since the scattering amplitude must still be dimensionless but now the amplitude involves a factor XjFf. the 
numerator must involve four powers of energy-momentum or two powers of energy-momentum together with rn^ or 
TO^. Thus any counterterm which is included in order to absorb this divergence must be /oitr-derivative in character. 
Gasser and Leutwyler studied this problem and wrote the most general form of such an 0(p‘^) counterterm in chiral 
SU(2) asi 


£4 = 


+ 


2=1 


tr(£)^ [/£>'" C/I') 




tr (xC/t + Ux^) + ^4tr [D^U^D^^x + D^UD^^x^ 


i/str (F^^C/F«'"^C/t) -t ^/gtr {Fff^D>^UD''U'^ + D''U) 




tr (x'I'C/ - Ux^) 
where the covariant derivative is defined via 


n 2 


( 22 ) 


D^U = d^U+ {A^,U}+[V^,Uf (23) 

the constants = 1,2,...7 are arbitrary (not determined from chiral symmetry) and are external field 

strength tensors defined via 


pL,R 


— 8 F^'^ — 8 


i[F. 


L,R jpLjRi 


pL,R 


= Vf,±A^. 


(24) 


Now just as in the case of QED, the bare parameters £i which appear in this Lagrangian are not physical quantities. 
Instead the experimentally relevant (renormalized) values of these parameters are obtained by appending to these 
bare values divergent one-loop contributions having the form 


(-i + 


7» 

327r2 


1 

e 


7 -I- ln( 47 r) -I- 1 


(25) 


where Xi are calculable constants^ and, evaluating the loop integrals in d dimensions, e = (d — 4)/2 and 7 is the 
Euler constant. By comparing with experiment, Gasser and Leutwyler were able to determine empirical values for 
the seven While seven sounds like a rather large number, this picture is actually quite predictive (T^. [Tlj| . 


3.2. Full effective field theory for QCD 

The use of effective Lagrangians is not the full content of an effective field theory. These do express some symmetry 
restrictions on amplitudes and also describe the residual effects from high energy. But effective field theory is a full 
quantum field theory with trees and loops of the particles that are described in the Lagrangian. The renormalization 
described above is only a part of the effects of loops, and not even an interesting part of the loops. Much more relevant 
for physics are the low energy effects of the loops. 

The actual calculations are straightforward and in the results we see several effects. There are the divergences 
described above, and these are absorbed into the renormalized parameters of the Lagrangian in a specified renom- 
rmalization scheme. Then there are these parameters themselves. Here we are faced with various possibilities. We 
can measure the parameters in other reactions, or perhaps turn to lattice calculations in order to predict them from 
the full theory of QCD. These numbers themselves are not predictions of the effective field theory. The distinction 
here is that the effective field theory has its structure determined by the symmetry and light degrees of freedom only, 
and can serve as the low energy limit of any theory which has these features. For example, the linear sigma model 
and also QCD-like theories with different numbers of colors all have the same structure for the effective field theory. 
However, the values of the constants in the Lagrangian would differ for such theories. In this sense, the constants 
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encode the physics of the ultimate high energy theory (often referred to as the “UV completion” of the EFT). Because 
these terms are local, this is what is expected from the uncertainty principle arguments discussed earlier. 

So what are the quantum predictions? We have argued that they are not the renormalization procedure nor the 
parameters themselves. However, there are residual quantum effects that are independent of the renormalization and 
of the parameters. These come from low energy propagation of the light particles. (The high energy parts of loops 
are local and respect the symmetry, so they contribute shifts to the renormalized parameters.) In advance we can 
know some things to look for. The local chiral lagrangian is expressed in powers of the derivatives. The resulting 
amplitudes in momentum space are then polynomials in the energies involved - these are an analytic expansion of the 
amplitude. However loops also bring in non-analytic terms, such as ln(—g^) or \/—q^, which cannot arise from an 
expansion of a local Lagrangian. These non- analytic terms are signals of long distance propagation by light particles. 
An important byproduct is that such non-analytic terms at one loop are always independent of any of the parameters 
of the Lagrangian as well as being finite and divergence free. 

In pure pion physics one finds uniquely the logarithmic non-analyticity, as this arises in the bubble diagrams that 
appear in the pionic theory. The square root \/—q^ arises in triangle diagrams with one massive and two massless 
particles, and then appears in pion interactions with baryons. In real QCD with massive pions, these non-analytic 
terms also appear as In(m^) or Such mass dependence is also a unique prediction of the effective field theory. 

The results in the third column of Table I are the results of this program carried out to two loop order @ also 
taking into account dispersion relations constraints, which is probably the gold standard for the pionic effective field 
theory of QCD. However there is an extensive literature of one loop studies of many processes. The full review of this 
program is not appropriate for the present document, but we should note that the program has been quite successful. 
We refer the reader to the literature for further information. 


4. EFFECTIVE FIELD THEORY OF GRAVITY 

Having seen how the strong interactions can be described via effective field theoretic methods, we move on to our 
primary goal, which is to treat the grav itational interaction in a parallel fashion, as developed in refs. [l^ . Reviews of 
this procedure can be found in [iSj, . The EFT method relies on an expansion in energy-momentum or derivatives 
and the gravitational interaction is described in this way. In order to understand how this is done we require the 
connection coefficient E^^, which is defined via the covariant derivative Da in terms of its operation on a vector field 
as 


DaA^ = daA^ + r^pA>^ 


(26) 


where, in terms of the metric tensor g^,^, 

ra/3 = [dagpc + dpgaC “ ■ (27) 

Taking the metric tensor as a field, we see that the connection involves a single field derivative, while the curvature 
tensor 


defined as 

\Dfi, Di,\Aa = Rafi^Ap 

has two, as do its associated quantities, the Ricci tensor 

D _ f) ■pQ! _ f) po; _j_ pC po: _ pC pen 


and the scalar curvature 


R = g^''R^,. 

The gravitational action can then be written as a derivative expansion 


S. 


grav 


= / 


A H— —R -\- ciR -l- C2RuvR^’^ -I- ... 


(28) 

(29) 

(30) 

(31) 

(32) 
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where A is the cosmological constant and k is given in terms of the Cavendish constant G via = 32ttG. The 
expansion parameter here is the Planck mass Mp = G~"s ^ 10^® GeV. Thus the higher order terms in the derivative 
expansion are suppressed by powers of E/Mp,p/Mp <« 1 so that the Einstein action should provide an extremely 
precise picture of gravitational effects at presently relevant energies[l^. Of course, a mystery is why the leading 
term in the expansion, the cosmological constant, is so small experimentally—A ^ GeV"^ [l^. (Note that a 

corresponding situation exists in QCD, where the smallness of the experimental result for the theta-term —9 ^ 10“^^— 
remains unexplained [ij.) However, we shall merely take this result as an empirical fact and will neglect A for the 
remainder of this paper. The full action is then given by including the matter Lagrangian Cmat- Varying the lowest 
order action 


Stot — Sa 


T Sjnat — 




Cri 


(33) 


we find the Einstein equation 


Rfj.!^ - 29^J.l'R = -SnGTfj,^ 

where 

T _ 

ay — / ^mat 

is the energy-momentum tensor of matter[l^. 

The theory can be quantized using the background field method by defining 


(34) 


(35) 


9tJ.u{x) = g^i.v{x) -I- Khf^u{x) 

where gfi^{x) is a classical solution of the Einstein equation. The inverse metric tensor is then given by 

= g^^'' - Kh^^^x) + + ... 


(36) 

(37) 


where indices are raised and lowered with the classical metric tensor . We choose to quantize about flat space, so 
that g^L, = rjfi,y. We find then for the Ricci tensor and scalar curvature at one derivative order 


Rp} = [dpd„h + 


dfi 


{d ■ h)i, - -d„h ] - { {d ■ h)f^ - -d^h 


1 


i?(i) = Kph-d^d^h^'^] 


Uh - 23'^ ( [d ■ h)a- - -dah 


(38) 


where h = {d ■ h)n = d^hx^, and □ = As is well known the Einstein equations are invariant under 

a general coordinate transformation which, in terms of the fields, implies a gauge invariance 


x^^ ^x^ = x^^+e^{x) 

^pv{x') t h ^iy(x) = duC^x') (39) 

for infinitesimal e^{x). In order to deal with this invariance, we must make a gauge choice and we elect to work in 
harmonic or deDonder gauge—= 0—which reads, to first order in the field expansion, 

0 = - ]^dah = {d ■ h)a - ^dah (40) 

whereby the linearized Einstein equation 




- ^gpuh) - df, ^(5 • h)„ - 
- ^{d ■ h)^ - + Vpvd^^ ^(5 • h)a - 


_ 2rpmat 

- 


(41) 









becomes 


□ I I = 


or its equivalent form 


r\h — _K- I _ ri 7^™“* 

\ ^nv 


(42) 


(43) 


where = 77^‘^T^“h There exist, of course, well known solutions to Eq. (j43|) . For example, in the case of a 

stationary point mass m located at the origin we have r^“*(a;) = 'q^io'qi,Qm5^{x), for which the solution of Eq. (H51) is 




with 


fi.r) = -h 


m 

167rr 


G m 
Stt r 


At next order it is useful to write the Einstein equation as 


where we identify 


= -f + T^/r) 




2 

as the energy-momentum tensor of the gravitational field. Using 


- ^dahfj,xd°'h^+ ^dahij,xd^h° 

+ [9\di,h^a + d\d^Ka - d^d^hxa - dxdahf,^,] 


+ 2 ~ ^ - dah^v) 


+ 2^ 


' 2dx ( {d ■ h)^ - ) - Uhxc. 


1 


-I- (9 • h)“ — -d°‘h ) (9 • h)ci — -dah 


1 


we find 


(44) 

(45) 

(46) 

(47) 


(48) 


= -‘^h^‘^{d^d,hx. + dxd,h^,-d,{d,h^x + d^Kx)) 

- 2dxK,d^hl + 2dxK,d'^hl - d^Kxd^h’^^ 

- Vt^iy {dxKxd^'h^^ - ~ - h^^ah 


Before proceeding we also need to deal with the gauge-invariance by using the Faddeev-Popov methodjl'i 
leading to a second order action of the form 


(49) 


= / d'^x 


-^df,haxd^h°‘^ + dah^^xd 


- h^“n/lAa - {d^'h^a - ]^dah 

H“ Sghosti'HfJ.) 


(50) 
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where is a fermion ghost field. 

Using the harmonic gauge condition Eq. (1401) and liberally integrating by parts we find 

^tot — 4 / ^ (x) ^(x) Sghostijlfi^ 


d^xh^’' {x)Pg„^al30h°^^ (x) + SghostiVfi) 


where 


Inverting, we see that the harmonic gauge graviton propagator is given by 

'^P^u,aj3 


dr gii>^a^{,d) — 


+ ie 


(51) 


(52) 


(53) 


With the propagator in hand we can begin to explore quantum gravitational effects provided we know the interaction. 
In the electromagnetic case we have the equation 


UAg + d^{d,A'^) = -eJg 

which, using the Lorentz gauge condition— d^A’^ = 0—corresponds to the interaction Lagrangian 


(54) 


CZt{x) = -eJg{x)A^{x) (55) 

Likewise in the case of the gravitational interaction we have the field equation Eq. dm which, using the harmonic 
gauge condition— d^hpg, — ^dgh = 0—corresponds to the interaction Lagrangian 


zrix ) = -T,.{xW%x) 


(56) 


so that the gravitational ’’charge” and ’’current” are k/ 2 and T^i, respectively. It is clear then that in order to 
determine the graviton couplings, we need to know the matrix elements of the energy-momentum tensor. 

We begin by considering the case of a scalar field, for which the energy-momentum tensor is derived in Appendix 
A—the lowest order scalar energy-momentum vertex is given by 


< P 2 \tIS){. x)\pi >= 


gi(P2-pi)-x 

'J AE 1 E 2 


2PgP, 




(57) 


where P = ^{pi + P 2 ) is the mean energy-momentum and q = pi — P 2 is the momentum transfer. Note that the 
energy-momentum tensor is conserved in that 


q^ < P2\Tf,y{x)\pi >= 0 (58) 

as required by taking the divergence of the field equation Eq. (1421) and using the harmonic gauge condition. Of course, 
radiative corrections—gravitational or electromagnetic—lead to modifications of the lowest order matrix element Eq. 
dSB, which must have the general form 

pi(p2—pi)-x 

< P2\T),p{x)\pi >= [2Pfj,P„Fi{q'^) + {q^q^ - ?7;xi.9^)7^2(g^)] (59) 

As can be seen from the condition 


< P2\Pii\Pl >= Pp < P2\Pl > = < P 2 I 


d^xTgo{x)\pi >= PpFi{q^ = 0) < p2bi > 


(60) 


conservation of energy-momentum requires 


Ei(g2=o) = l (61) 

but there exists no constraint on We see from Eq. (I57p that at lowest order Ei = 1 and E 2 = — ^, so that this 

condition is satisfied. 
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5. PHOTONIC LOOPS FOR THE REISSNER-NORDSTROM AND KERR-NEWMAN METRICS 


Before examining the modifications induced by gravitational corrections, it is useful to examine a related but simpler 
case—electromagnetic radiative corrections to —in order to understand the relevant physics [l^ [l^. Calculation- 
ally the photon loop corrections are very similar to those arising from graviton loops because photons and gravitons 
are both massless and propagate long distances. The example of the photon will introduce the nonanalytic corrections 
that occur in momentum space and will show how these are correspondingly nonlocal in position space. 

Our subject is the photon loop correction to the energy-momentum tensor of a charged particle, which will reveal 
the form of the gravitational field in the vicinity. The classical result should be the Reissner-Nordstrom metric and 
Kerr-Newman metric, which describe the gravitational field around a charged particle without and with spin. Because 
of the charge, the photon loop is the leading correction. Gravity can be treated classically here. 

In treating the photon loop, we need the gravitational coupling to the photon. In lowest order this involves the 
electromagnetic energy-momentum tensor 

(62) 

which leads to a lowest order energy-momentum tensor for the photon 


{p2,e2\Tfj,,,{x) \pi,€i) 


-f 


-h 


ai(p2-Pl)-X 


Pp (eLei • 1 ' q) + Pv ’ 9 “ ' «) 


2 ei • 62 - -qe^-q 

* , Q.v ^ * 

Ye 2 ^ei • q + ■ q + ' 9 + '9 

y / I \ 

~2 


(63) 


Eq. (|63)l defines the Feynman rule for the photon coupling. 

We can now evaluate the various photon loop diagrams, which will lead to modifications of the lowest order spin 
zero form factors. The scalar electromagnetic vertices are well known and, using the photon-photon-graviton vertex 
given in Eq. (l63)) . the results arefl^ 


0 ^em(^ 2 ) 


3 TT^m 



91 ^ 

2 log - - log - 

771.^ 3 m 



47r 


-- -f 7- 
e 


log 


47r^2 



(64) 


where (Xem = e^/ 47 r is the fine structure constant. Here we have used dimensional regularization and A is a pho¬ 
ton “mass” which is inserted in order to regulate the infrared sector of the theory and which disappears when 
bremsstrahlung corrections are included. We observe that there exist two types of radiative corrections here. One 
class is analytic and therefore local when the transition to coordinate space is made. (This includes the ultraviolet 
divergence, which can be absorbed into the coefficient of a term RF^^F^'^FiQU in the effective Lagrangian.) 
The second class is more interesting and involves nonanalytic terms such as \/—q^ and q^ log —Such forms do not 
occur in radiative corrections to the electromagnetic current and arise here from the triangle and bubble diagrams in 
Figures la and lb involving coupling of the energy-momentum tensor to a pair of photons. It is the presence of the 
two massless propagators in such diagrams, and the fact that both photons can be nearly on shell, which leads to this 
nonanalytic structure[ 2 l|. 

The physics of these nonanalytic terms can be easily extracted by making the transition to coordinate space in the 
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(b) 


FIG. 1: Electromagnetic correction diagrams having nonanalytic components. Here the single wiggly lines represent photons 
while the double wiggly line indicates coupling to a graviton.. 


Breit frame, wherein go = 0 and pi = —p 2 = q/2— 

cPq 


0 


or-(r) = 




(27r)3 

cPq 


+ F2{-q^) 


^iqr 


TTQle 


m — 


aemq^ . 

^-log- 

OTrm 777 


Orpem/ \ _ 

^ ’ 2m 
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0 

1 /■ (Pq 


+ ... 




ttem If \ Oem^ 

dTrr"^ I 2 


f TTCTem tre772 g 

\ 16|q| Sirm ^ m? 

+ ... 


(65) 


where we have restored the factor of h in order to differentiate classical and quantum mechanical contributions and 
the ellipses denote short distance pieces. The interesting feature here is that a quantum loop diagram has generated 
a classical effect—that is, a term independent of h. As mentioned above this is due to the presence of two massless 
propagators in the diagrams [2lj. The meaning of these classical pieces is quite clear and can be understood by using the 
classical electromagnetic energy-momentum tensor which results from the lowest order electromagnetic Lagrangian. 
In the vicinity of a stationary particle having charge e and located at the origin we have E{r) = ef/47rr^ so that 


clT^rn^r) 


2 STrr^ 

0 

—EiEj -\- -SijE"^ = 



( 66 ) 


which agree precisely with the classical components of Eq. (1651) . The meaning of the /i-dependent corrections can be 
understood qualitatively by including the effects of zitterbewegung. The point is that at a classical level the distance 
of the probe from the particle can taken to be a fixed distance r. However, including quantum mechanical effects, 
the location of the source is uncertain by an amount of order the Compton wavelength— 5r ^ h/m « r. Thus l/r^ 
should be replaced by the form l/(r -|- 6r)'^ ^ which has the form of the quantum corrections found above. 

The feature that the loop correction leads to well-understood classical corrections is also valid if the source particle 
has spin. In the case of spin 1/2 the form of the radiative corrections to the energy-momentum tensor have also been 
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calculated. In this case, as shown in Appendix A, the lowest order energy-momentum tensor vertex has the form 
< {x)\pi >= - to)] u{pi) 

After radiative corrections the energy-momentum tensor can be written in the general form 


(67) 


< P2|7)ji.(a:)|pi > = ^^^''^u{p2) 


—P^P^ ^^Fi{q^) + —{q^^qu - = ^ 2 ( 17 ^) 

m m 

i 


so we see that in lowest order we have 


u{pi) 


1^(0) ^(0) ^ 1 1^(0) ^0 


( 68 ) 


(69) 


In the spin 1/2 case, besides the constraint of energy-momentum conservation discussed above we have the additional 
requirement of angular momentum conservation. Defining 

Mi2 = J d^x{ToiX2 - T 02 X 1 ) -iVg, j dPxP'f^Toii.r) + J d?xPf^To2{r), (70) 


Since 


^ = 1 ™ <7'2,t 1^12^1, t>= Mt(p)^CT3Ut(p) = F3(g^ = 0) 


(71) 


we see that angular momentum conservation requires that ^ Fs^q'^ = 0) = 1. Again energy-momentum conservation 
requires 2 Fi{q^ = 0) = 1 while there is no constraint on 2 ^ 2 ( 9 ^). Obviously the lowest order forms given in Eq. (IMl) 
satisfy these conditions. 

Performing now the loop integrations associated with electromagnetic corrections, we determine that 


^E/™(g 2 ) =0 F/™(g 2 ) and ^F^'^iq^) =° F^^iq"^) 
while for the new form factor iF§™{q'^) 


(72) 


^Fi'^{q^) = 1-f 
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For the spin 1/2 energy-momentum tensor we find 
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(74) 
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where S = That is, the diagonal spin zero and spin 1/2 energy-momentum tensor densities are identical— 


iTSS^ir) =° TSS^ir) and =° 7/7(r) 

—but there now exists a nonzero off-diagonal term The form of the classical correction in this term can be 

understood from the feature that a Dirac particle has both an electric field and magnetic field 


E = 


ef 

47r7’2 


and 


e 3fS ■ f — S 
m ‘iirr^ 


(75) 


The off-diagonal piece of the classical energy-momentum tensor density, Eq. (ESI) 


= -{E X B), 


(76) 


then becomes 




47rmr® 


{S X r)i 


(77) 


in agreement with the classical component found in Eq. dZll). The form of the corrected energy-momentum tensor 
has also been calculated for a spin 1 particle, yielding identical results as for the spin 1/2 case, except for the 
replacements XfX^ ^/ ’ ief x ii plus new quadrupole corrections [^. We suspect that the forms 

of these corrections—both classical and quantum—are universal. 

Before moving the gravitational case it is useful to note one other interesting feature in the electromagnetic case. 
Using the linearized Einstein equation—Eq. 03]— and the results for the energy-momentum tensor T^^ir) generated 
above we can solve for the metric tensor, which yields the form 


°hoo{r) 

°h,o{r) 


for a spinless particle and 
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(78) 


(79) 


for a particle with spin 1/2. (The results for spin 1 have also been calculated and agree with those of spin 1/2 up to 
small quadrupole corrections [2^ so again, it is likely these results too are universal]^.) The classical components of 
the spin zero results Eq. (175)) agree precisely with those of the Reissner-Nordstrom metric [2^, which is the metric 
associated with a massive charged particle, while the spin 1/2 results Eq. (17^ agree with those of the Kerr-Newman 
metric (^. which is the metric associated with a massive charged particle which is spinning. Again the form of the 
quantum mechanical corrections are consistent with zitterbewegung fluctuations. 


6. GRAVITATIONAL CORRECTION TO THE SCHWARZSCHILD METRIC 

We now repeat the same procedure but with gravitational loops. This amounts to looking at graviton loop correc¬ 
tions to the Schwarzschild metric. This procedure is not gauge invariant, and even the concept of a metric is not a 
fully quantum concept. However, the result is an illustration of the form of quantum corrections. We are working 
in harmonic gauge and the result applies only in that gauge. However, in the process of calculating the quantum 
correction, we also obtain the first classical correction to Schwarzschild, a result first found by Duff|^. 

Again there exist nonanalytic forms arising from the triangle and bubble diagrams containing two massless 
propagators—Figure 2—which lead to classical and quantum mechanical corrections to the lowest order results. Us¬ 
ing the gravitational couplings given in Appendix A and keeping only the the nonanalytic pieces, the gravitationally 
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corrected form factors in the case of a spinless particle of mass m are found to be 


» frv ) = i + 


TT \ 16 y/—q^ 




= -t + 2 ^U 4 ^- 21 og^ 


,-q^ 


(80) 


The corresponding energy-momentum tensor is 




= m5^{r) — 


3Gm? 3Gmh 
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• ry-M = 
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47rr^ 


r 2 2 ^*' 


2Gmh ^ 


(81) 


In this case we can compare the classical results with the predictions of the gravitational energy-momentum tensor 
density—Eq. (|4^ — 

3Gm^ 


'Tr^'ir) = (-3V/(r) • V/(r) - 12/(r) V2/(r)) = + .. 


Tfo™"(r) = 0 


Tfrir) = (-2Vg/(r)V,/(r) + V/(r) • V/(r) - 4/(r)VgV,/(r) 

T 4.g,/(r)VV(.))^-^(^4^.)T... 


(82) 


where the ellipses denote short distance components, which agree precisely with the classical component of the loop 
calculation result—Eq. m- 



FIG. 2: Feynman diagrams having nonanalytic components. Here the doubly wiggly lines represent gravitons. 


In the case of spin 1/2, the gravitationally corrected form factors are found to be 

Lr^grav, 2\ i ^ 3 , 
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That is, the spin zero and spin 1/2 values for the form factors Fi, Fj are identical— 


—which in turn implies that the diagonal components of the energy-momentum tensor densities are are identical to 
those of spin zero— 

°Foo(r) Foo(t') and °ry (r) =5 F^(r) 


However, there is now a nonzero off-diagonal piece 

= liS X V),S^{r) 



l5Gh \ 

47r2j^7 j 


{S X r)i 


(84) 


where S = which can be compared to the predictions from the classical gravitational energy-momentum 

tensor density—Eq. (|4^ — 




[-{S X V),/(r)V,V,/(r) + V,/(r)V,(S x V),/(r)] 


Gm 

27rr® 


{S X r), 


(85) 


Again there is complete agreement between the two classical calculations. (The spin 1 calculation has also been 
performed and agrees with the spin 1/2 forms, except for small quadrupole corrections [l^, so that we suspect that 
these results are universal.) 

As in the electromagnetic case, we can use the linearized Einstein equation to convert the results for the energy- 
momentum tensor to those for the metric tensor, yielding for spin zero 
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2Gm 2G^rn? IG^mh 
- 1 ^ 1 5 — 


= 0 
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( 86 ) 


the classical components of which agree completely with the Schwarzschild metric [2^, which characterizes a stationary 
massive particle. 

In the case of spin 1/2 the diagonal components of the metric tensor are identical to the spin zero case —2 /iq™’' (r) =° 
/iQQ“’'(r) and 2 hF“’'(r) =° —but there exists now a nonzero off-diagonal component. 



2G^m 3G^h\ 

-4-'-5“ ) 

G Trr° J 


{S X r)i 


(87) 


The classical components of the spin 1/2 metric tensor agree completely with the Kerr metric [l^, which characterizes 
a massive spinning particle. Again, the spin 1 result has also been calculated and agrees completely with the spin 1/2 
forms [ 2 ^, so we suspect that they are universal. 


7. CORRECTION TO THE NEWTON POTENTIAL 

It has long been thought to be an unattainable goal to calculate effect of quantum physics on the gravitational 
interaction. However, for the long distance quantum correction to the Newtonian potential, the result is both simple 
and universal. We explain the logic in this section. 

There are several possible definitions of a gravitational potential. We shall discuss some of the associated subtleties 
below, but for the moment we shall simply define the potential V (r) as the Fourier transform of the nonrelativistic 
scattering amplitude M{q )— 

and we begin our discussion by considering the scattering of a pair of spinless particles. 


(88) 
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a i, 

FIG. 3: Basic kinematics of gravitational scattering. 


7 . 1 . Spin 0-Spin 0 Scattering 


At lowest order the interaction of two spinless particles of mass mi, m 2 is described in terms of the one graviton 
exchange (tree) amplitude 


= -i 
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(89) 


We shall utilize the symmetric center of mass frame with incoming momenta pi = —p 2 = P E q/i and outgoing 
momenta pa = —P 4 = p — q/2. Conservation of energy then requires p - q = 0 so that p^ = p ^ -t- q ^/4 for f = 1, 2, 3,4 
and q^ = —q^. In the nonrelativistic limit— q^,p^ <C m^ —the lowest order amplitude reads 
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yielding the potential 


V«(r) = - 
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Gmim2 
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(91) 


and we recognize the Newtonian potential as the dominant piece of Eq. (1911) (accompanied by a small kinematic 
correction) together with a short range modification. 

Our goal is to examine corrections to this lowest order potential due to two-graviton exchange and thereby to 
define a higher order gravitational potential. This problem has been previously studied by Iwasaki using noncovariant 
perturbation theorvpq. and by Khriplovich and Kirilin[^.[30l| and by Bjerrum-Bohr, Donoghue, and Holstein}^ 
using conventional Feynman diagrams. Our approach will be similar to that used in and (3Tj |. The diagrams 

utilized are shown in Figure 3 and the various interaction vertices are derived in Appendix A so it is merely a matter 
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of calculating these triangle and bubble diagrams. Because of the many tensor indices involved, this is a challenging 
but straightforward problem and until recently there had been a number of mistakes in such evaluations . which 
have finally been corrected[^. As before, the procedure is to calculate the various diagrams while retaining only 
the nonanalytic components since only these terms lead to long range corrections to the Newtonian potential. The 
nonanalytic contributions which arise from the various diagrams can be expressed in terms of the quantities L = 
and S = 7 r^/-\/|q| Summing all the scattering diagrams, we determine the total 
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6 (mi + 7712 ) 6 '- —L 
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777 1 7772 
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(92) 


and observe that, in addition to the expected terms involving L and 6 , there arises a component of the second order 
amplitude which is imaginary and represents a scattering phase. The origin of this imaginary term is from the second 
Born approximation to the Newtonian potential, and suggests that in order to define a proper correction to the first 
order Newtonian potential we must subtract off such pieces. For this purpose we work in the nonrelativistic limit and 
the center of mass frame —pi + P 2 = 0—as dehned above. We have then 


s - So = 



27771 7772 


(93) 


and 


77717772 777 ^ 

S - So ~ Po 


(94) 


where nir = 77717772/(7771 + 7772 ) is the reduced mass and po = |Pi|) * = 1)2, 3,4. The transition amplitude Eq. (IMl) 
then assumes the form 
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We can evaluate the Born iteration directly by utilizing the simple Newtonian potential 


(95) 


(96) 


which reproduces the long distance behavior of the lowest order amplitude for spin -0 - spin -0 gravitational scattering— 
Eq. (I91[) — in the nonrelativistic limit. The corresponding momentum space representation is 


and the second Born term becomes 
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which precisely reproduces the imaginary component of °A4joj(q), as expected. In order to produce a properly defined 
second order potential ^Vq^v) we must subtract this second order Born term from the second order transition 
amplitude, yielding a well-defined second order gravitational potential 
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The quantum mechanical—^ h/mr ^—component of the second order potential given in Eq. (1991) a grees with that 
previously given by Bjerrum-Bohr, Donoghue, and Holsteinand by Kirilin and Khrir)lovich[29j. However, the 
classical—^ 1 /r^—contribution quoted by Iwasaki 


0 



G^mim2{mi + m 2 ) 


( 100 ) 


differs in both sign and magnitude from that quoted above in Eq. (|M)) and by Bjerrum-Bohr et al. in [^. The 
resolution of this issue has been given by Sucher, who pointed out that the form of the classical interaction depends 
upon the precise definition of the first order potential used in the iterationf^. Moreover, it depends and on whether 
one uses relativistic forms of the leading order potentials and the nonrelativistic propagator (£) in the iteration. 
In modern terms, the potential depends on how one performs the matching— e.g., Iwasakiperforms an off-shell 
matching while we match on-shell.^. Use of the simple lowest order form Eq. dSH) within a nonrelativistic iteration 
yields our result for the amplitude given in Eq. (j98p and is sufficient to remove the offending imaginary component of 
the scattering amplitude. In Appendix B we derive an alternative form of the G(G^) classical potential which results 
from an iteration that includes the leading relativistic corrections and which reproduces the Iwasaki result [^. 

Therefore, a unique definition of the second order potential potential does not exist. However, ambiguities in the 
form of the second order classical potential are not a concern, since the potential is not an observable. What is an 
observable is the on-shell transition amplitude, which is uniquely defined in each case as 


°Mtotiq) 
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( 101 ) 


where the index i denotes differing possible definitions of the potentials and the iteration. Thus we regard the potential 
as merely a way to display the resulting scattering amplitude in coordinate space, and we emphasize that the main 
results are the long distance components of the scattering amplitude—With these caveats in mind, the 
total potential describing the gravitational scattering of spinless particles, at second order in G, can be written as 
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( 102 ) 


and we observe that there exist long range contributions to the leading Newtonian potential, with a classical component 
falling as 1/r^ together with a quantum mechanical corrections dropping as h/mr^. 


7.2. Spin 0-Spin | Scattering 


The calculation of the correction to the Newtonian potential can also be carried out straightforwardly in the case of 
a spin 1/2 particle having mass m 2 scattering from a spinless particle of mass mi. The tree level transition amplitude 
from one-graviton exchange is 
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(103) 


Defining the spin vector as 

S 2 = -^uiPi)l'"'y5u{P2) ^ xlf^o■X2^ (104) 

the nonanalytic part of the transition amplitude in the threshold limit s —>■ sq = (wi + ^ 7 ^ 2 )^ can be written in the 
form 
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1 


Besides the dependence on the forms used in the iteration, the classical piece also depends on the coordinates used. The quantum piece 
however depends neither on the choice of coordinates isill nor on the iteration forms O . 
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In order to define the potential we again take the nonrelativistic amplitude in the symmetric center of mass frame 
(Pi = -P 2 =P + <7/2)— 
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and the lowest order potential becomes 
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(107) 


where L = rxpis the angular momentum—the modification of the leading spin-independent potential has a spin-orbit 
character. 

A subtlety that arises in the calculation involving spin is that two independent kinematic variables arise: the 
momentum transfer and s — Sq, which is to leading order proportional to Pq (where i = 1, 2, 3,4) in the 

center of mass frame. We find that our results differ if we perform an expansion first in s — sq and then in or vice 
versa. This ordering issue occurs only for the box diagram, diagram (d) of Fig. 2, where it stems from the reduction 
of vector and tensor box integrals. Their reduction in terms of scalar integrals involves the inversion of a matrix whose 
Gram determinant vanishes in the nonrelativistic threshold limit q^,s — sq —>■ 0. More precisely, the denominators 
or the vector and tensor box integrals (see Appendix A in [s^) involve a factor of (4pg — q^) when expanded in the 
nonrelativistic limit. Since q^ = 4pgsin^ ^ with 9 the scattering angle, we notice that dpg > q^ unless we consider 
backward scattering where 9 = tt and where the scattering amplitude diverges. And since pg originates from the 
relativistic structure s — sq, it is clear that one must first expand our vector and tensor box integrals in q^ and then 
in s — Sq. 

Calculating the various diagrams as before we find the total second order contribution 
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Finally, working in the center of mass frame and taking the nonrelativistic limit, Eq. (nns) becomes 
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We observe from Eq. (11091) that the scattering amplitude consists of two pieces— 
i) a spin-independent component proportional to X 2 /X 2 i whose functional form 
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is identical to that of spinless scattering. 
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ii) a spin-orbit component proportional to 
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whose functional form is 
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We note in Eq. (11111) the presence in the spin-orbit potential of an imaginary final state rescattering term proportional 
to i/po, similar to that found in the case of spin-independent scattering, together with a completely new type of 
kinematic form, proportional to l/p§ which diverges at threshold. The presence of either term would prevent us from 
writing down a well defined second order potential. 

The solution to this problem is, as before, to properly subtract the iterated hrst order potential— 
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We find that the iterated amplitude splits into spin-independent and spin-dependent pieces. The leading spin- 
independent amplitude is 
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and the leading spin-dependent term is 
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so that when the amplitudes Eqs. (I116p and (11151) are subtracted from the full one loop scattering amplitude Eq. 
(I109|) both the terms involving 1 /pq and those proportional to i/po disappear, leaving behind a well-defined second 
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We observe that the second order potential for long range gravitational scattering of a spinless and spin-1/2 particle 
consists of two components: one independent of the spin of particle 2 and identical to the potential found for the case 
of spinless scattering, accompanied by a spin-orbit interaction involving a new shorter range form for its classical and 
quantum components. 


7.3. Universality 

We have seen that the leading quantum correction to the Newtonian potential is independent of the type of particle 
being considered. This result obtains despite the fact that very different Feynman diagrams occur when dealing with 
fermions and bosons. In fact, one can prove that that this universality is itself a low energy theorem of quantum 
gravity [ 3 ^. 

The argument which demonstrates that this equality is more than an accident has been shown by the use of some of 
the new methods of quantum field theory, which also serve as a check on the Feynman diagram calculation. The new 
techniques are often referred to as unitarity methods [s^, because they rely on the unitarity and analyticity properties 
of Feynman diagrams. They work by identifying the unitarity cut in the amplitude and reconstructing the full 
Feynman diagram from this information. All one loop Feynman diagrams can be reduced to scalar bubble, triangle 
and box diagrams without any factors in the numerator of the loop integral, a property referred to as Passarino- 
Veltman reduction. Each of these structures has distinctive unitarity cuts. From the cuts then, one can reconstruct 
the prefactors of the box, triangle and bubble diagrams. In addition there can be polynomial terms which do not lead 
to cuts but, as we have above, we are interested only in the nonanalytic terms, which can be reconstructed properly. 

For the calculations described above, we need only take the gravitational Compton amplitude—the coupling of 
two on-shell gravitons to two on-shell matter particle—and multiply them together in order to get the two-graviton 
unitarity cut. The contraction is performed most simply using helicity methods [3^. which involve a form of axial 
gauge. 

A second modern miracle further simplifies the amplitude method. It has recently been discovered that on-shell 
gravity amplitudes are in a precisely specified way related to on-shell gauge theory amplitudes [s^. This property is 
summarized by the phrase: gravity is the square of a gauge theory. In our cases, the gravitational Compton amplitude 
is the square, with a given prefactor, of the QED Compton amplitudes. The gravitational amplitude is very complex 
because of the presence of the triple graviton vertex, but the QED analog can be worked out straightforwardly by 
any field theory student. 

These calculations have been carried out[35l|. reproducing the results of the Feynman diagram approach. This is 
gratifying as it shows both that the calculations have been done correctly and also that they are gauge invariant, as 
the two methods use different gauges. However, they also provide a proof in the universality low energy theorem. This 
is because the Compton amplitudes are already known to have universal soft limits [3^. In the product, the leading 
term is then universal - this is the one that gives the classical correction. The quantum result follows from a term 
in the product which is a factor of \/—q^ higher than the leading term. However, the non-universal features appear 
only at a power q^ higher than the leading term. This implies that the leading quantum correction is also universal. 

It is worth mentioning that both the classical and quantum corrections have also been reproduced by dispersion 
relations methods, as reported in [35j |. This is another technique which reconstructs the real parts of the amplitudes 
from their on-shell cuts. The non-analytic terms that we are interested in are independent of the number of subtractions 
needed. The dispersive method has been done in both the harmonic gauge, in which case one needs to include cuts 
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from ghost fields, and in the axial gauge of the helidty method, which has no ghosts. Again, this calculation can be 
used to prove the universality low energy theorem. 


8. GRAVITATIONAL SCATTERING OF A MASSLESS SYSTEM FROM A MASSIVE SYSTEM 

The calculations described above have dealt with the influence of gravity on massive systems and the results were 
based on an expansion in powers of momentum transfer over mass. However, gravity also couples to massless systems 
such as the photon. Here the calculations are more complicated because both the photon and the graviton can 
propagate long distances in loops. The first published calculations on these systems were described in and used 
the unitarity based methods described in the previous section. However, there is also an unpublished thesis [4l| 
which had studied the same systems^ using conventional Feynman diagrams. We will use the notation from both 
descriptions in what follows. 


8.1. Interactions of a massless scalar 
For simplicity we begin with the case of a (fictitious) massless scalar field. 

Consider a massless scalar particle of energy E which is moving along the z-axis. The lowest-order energy- 
momentum tensor of a massless scalar particle is still given by Eq. (EZl, but now with Ex = |pi|, E 2 = \P 2 \- 
For the particle to act like a fixed source, pi ~ p 2 ~ P >> Q, and so 

{P 2 \T^^'^ (x) \pi) ft! (118) 

where E ft Ei ft if 2 is the time component of P = ^(pi -I- P 2 ) and Further, 
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SO 

q^=Ex-E2-^-q = q., ( 120 ) 

IpI 

and q ■ X ft qz (t — z) — q^x — qyp, where the z-axis is taken to be in the direction of P. Then, transforming to 
coordinate space, 

(x) = [ =E6{z-t)d (x) S (y) (121) 

J (27r) 

which is the Aichelburg-Sexl form for the energy-momentum tensor of a massless particle of energy 


8.2. Metric Tensor 


The corresponding metric tensor is given by solving the (linearized) Einstein equation 

Ohyy{q) = -IQttG < P2\Tfj.u{x)\pi > 
where we have used the fact that TrT = 0. We have then 
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^ A few mistakes in the thesis were corrected by f^ . 
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which is the form of the metric given by Aichelberg and Sexljd^. 

In order to calculate the gravitational loop corrections to Eq. (11211) . we must evaluate the same diagrams as in the 
massive case and the result can be written in terms of the two form factors Fi{q^)^ F 2 {q^) defined in Eq. (l59l) . yielding 
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where the particle mass has been used as a regulator. Notice that there exist no classical nonanalyticities— 
^ \/—q^ —here. This result is obvious in retrospect since there exists no mass scale to divide by and is consistent 
with the feature that Aichelberg and Sexl demonstrated explicitly that the linearized solution, Eq. (I123|) . is also a 
solution of the full Einstein equation—there exist no higher order classical contributions . 

To lowest order in q (and therefore at longest range) the energy-momentum tensor in co-ordinate space can naively 
be obtained by taking the Fourier transform, yielding 
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However, there exists an obvious problem here in that Eq. (11251) depends on the artificial graviton mass A. Unlike the 
case of a massive particle, where the A-dependence appears only in a short range term, the A-dependence in Eq. (|125D 
is contained in a long range component and therefore we conclude that the long range energy-momentum tensor is 
not an observable. This A-dependence is not a problem, however, for the scattering amplitude, which is an observable 
and is the quantity which we next examine. 


8.3. Massless Particle Gravitational Scattering 

Consider now the gravitational scattering of a test particle of mass m (later taken to be massless) from a heavy target 
mass M, both of which are taken to be spinless. The needed diagrams are identical to those required for evaluation 
of the massive case, with the exception that we now must include the contributions from the bremsstrahlung terms— 
both from the massless (m) and massive (M) particles—which were unimportant in the massive scattering case since 
they were associated with short distance (analytic) effects. However, there is an additional feature which must be 
addressed in the massless scattering situation, which is the prevalence of infrared singularities. We know from QED 
that there are soft singularities which arise when loop momenta get small. Gravity has these soft singularities also. 
In Yang-Mills theories with massless charged particle, there are also “collinear divergences” , which lead to factors 
of log m? and arise when one of the external momenta of the massless particles is parallel to a loop momentum. In 
QED this collinear effect could in principle arise but, since there exist no massless charged particles, such m —>■ 0 
divergences are not an issue. In gravity, one might think that these singularities also exist as the massless gravitons 
carry their charge. However Weinberg [44l| showed that gravity does not have collinear singularities, and so all that 
we need to deal with are the soft infrared divergences. 


8.4. Massless Particle-Massive Scalar Scattering: Result 


We begin with the gravitational interaction of a heavy scalar of mass M with a light scalar of mass m, which will 
be later taken to vanish. The elastic differential scattering cross section is then given by 


dcTe; 


(27r) S {pi +P 2 -P 3 + Vi) 


\Mt 


d Pad p4 


4y^(pi -^2)^ - MW 2E3 ( 27 r)^ 2A4 


(126) 


The calculation must be done very carefully and details are given in However, the use of traditional Feynman 
diagram methods is quite tedious as well as challenging, and there exist a few errors in this result. As an alternative 
then, Bjerrum-Bohr et al. performed the evaluation using modern on-shell helicity amplitude techniques (4^, wherein 
the calculation is greatly simplified by at least three features: 
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FIG. 4: The two graviton cut for the amplitude between a massless particle (dashed line) and a massive scalar (solid line). The 
grey blob are tree-level gravitational Compton amplitudes. 


a) One is that the on-shell gravitational tree-level amplitudes can be written as the square of gauge theory am¬ 
plitudes [ 4 ^ . In the case at hand the (nonabelia n) g ravitational Compton amplitudes are reduced to the 

product of (abelian) QED Compton amplitudes [ssl l47l 1^ . The challenging diagrams involving the triple 
graviton vertex are avoided and are replaced by much simpler evaluations involving only QED vertices. The 
general relation connecting the gravitational and electrodynamic Compton processes is derived in detail in (48l | 
and is given by 
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where h represents a graviton while rj can be either a photon 7 or massless spinless particle ip. Here = 

■^ 5 = 1 * is th® Compton amplitude for the scattering of a photon from a massless charged spin-1 target while 
_^QED _ ^QED j.gpj.gggj^^g ^j^g Compton amplitude of a photon from a massless charged spin-0 target. These 
tree-level relations connect one-loop gravitational physics with one-loop electrodynamics in a non-trivial and 
interesting way [s^ . 

b) The second great simplification involves the use of on-shell unitarity techniques [i^, instead of Eeynman dia¬ 
grams. Unitarity-based calculations construct the relevant amplitude from the discontinuities of the scattering 
amplitude. The long range nonanalytic terms in the one-loop amplitude can then be readily calculated from 
these on-shell cuts using the property of unitarity, as was directly demonstrated in ref. [ 3 ^ . Cutting the graviton 
internal lines, the integrand of the one-loop amplitude factorizes in terms of a product of relatively simple tree 
amplitudes, given in this case by the gravitational Compton amplitudes. 

c) The final simplification is the use of the spinor-helicity formalism (see for a review). While this notation is 
perhaps less familiar, it drastically simplifies the form of the amplitudes which we display. 


The tree-level massive scalar-graviton Compton amplitude is 
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The tree amplitudes connecting a massless scalar ip and the graviton are then obtained by taking the limit M —>■ 0. 
Amplitudes with opposite helicity configurations are obtained by complex conjugation. Computation of the cut 
discontinuity can be accomplished using traditional methods and is greatly simplified by the use of the on-shell 
identities. 
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Using these techniques the leading contribution to the amplitude (expanding all integrals in terms of leading order 
contributions as done in [H, is found to be: 
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where /i^ is the arbitrary mass scale parameter used in dimensional regularization and hu^ = 3/40. 

In Eq. I129l the two terms in the second line correspond, respectively, to the leading Newtonian contribution and first 
post-Newtonian correction [^. [^. The next three (logarithmic) terms represent quantum gravity modifications. 
The first term on the third line corresponds to the quantum correction to the metric evaluated in (53j |. The second 
piece on the third line arises from the one-loop ultraviolet divergence of the amplitude and is the only contribution 
depending on the spin of the massless field. On the fourth line the first term involves a new form not found in the 
previous (massive) analysis. Finally, the last term, arising from the discontinuity of the box integral, contributes to 
the phase of the amplitude and is not directly observable. For this reason it will not be considered further. 

One can straightforwardly generalize the calculation to the case that the massless scalars are replaced by photons. 
In this case the only nonvanishing gravitational Compton helicity amplitudes involving photons 7 and gravitons h are 
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with by the above formula with pi and p 2 interchanged, and amplitudes with opposite helicity 


configurations are obtained by complex conjugation, 
massive scalar and photon is then found to be 


The resulting gravitational Compton amplitude involving a 
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where bu^ = —161/120 and = ( 2 Ma;)^/( 2 (pi|p 3 |p 2 ]") for the (- 1 —) photon helicity contribution and its complex 
conjugate for the (—h) photon helicity contribution. The photon amplitude vanishes for the polarization conhgurations 

(-1—1-) and (-), which is a direct consequence of the properties of the tree-amplitudes in eq. (I130L 

In contrast with the non-relativistic case, where there was a universality theorem for the coefficient of the quantum 
correction, we see that most terms agree, except for one. Comparing Eqns. 11291 and 11311 the exception is seen to 
be the bu^ contribution from the massless bubble diagram (77 = We note also that, because of the 

vanishing of the photon scattering amplitudes for the helicity configurations (- 1 —h) and (-), the amplitudes parallel 

and perpendicular to the plane of scattering are identical, which rules out the existence of birefringent effects. 


8.5. Significance 

There are some quick and general conclusions that can be drawn from this calculation. One is that massless particles 
no longer move along null geodesics. We have obtained the classical behavior as a first approximation, and indeed 
we even have obtained the correct result to second order. However there are new effects the interaction which have a 
different power dependence which will then modify the trajectory. Moreover, some of these interactions are different 
depending on the type of massless particle, a scalar vs a photon. This is a violation of some classical forms of the 
equivalence principle. The equivalence principle itself means different things in different settingsjs^, yet classically 
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would state that massless particles move on null geodesics which are the same for all massless particles. Both the 
geodesic motion and the universality turn out to be violated by quantum corrections. Both of these effects can be 
ascribed to the long distance propagation of particles in loops. These loop diagrams then sample spacetime points 
different from the classical geodesic. Particles have intrinsic power-law nonlocality and this leads to non-classical 
motion. 

The most familiar application of the scattering of massless and massive systems is probably the bending of starlight 
by the sun. A fully quantum treatment of this light bending which is capable of including the one-loop amplitude 
effects is not available. However, in order to try to understand the impact of these corrections, one can proceed by 
defining, in the small momentum transfer limit t ~ — 9 ^, a semi-classical potential for a massless scalar and photon 
interacting with a massive scalar object by use of the Born approximation result 
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where rg is an infrared scale. 

Using naively the semi-classical formula for angular deflection given in 
we find the bending angle of a photon and for a massless scalar 
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63 


The first two terms give the well known classical values, including the first post-Newtonian correction, expressed in 
term of the gauge-invariant impact parameter b (see for instance 57|). The last term is a quantum gravity effect 


of order G^hM/P = i'j^rs/{2b^) and involves the product of the Planck length and the Schwarzschild radius of the 
massive object divided by the cube of the impact parameter and depends on the spin of massless particle scattering 
on the massive target. Of course, this dependence does not necessarily violate the equivalence principle in the most 
fundamental sense, in that the logarithmic quantum corrections correspond to non-local effects in coordinate space. 
Because of quantum loop effects, the long-distance propagation of massless photons and gravitons is not localized, 
and consequently can be interpreted as a tidal correction in that the massless particle is no longer be describable as a 
point source. There is then no requirement from the equivalence principle that such non-local effects be independent 
of the spin of the massless particle. 

Numerically, we can compare the bending angle of a photon with that of a massless scalar by the sun. The only 
difference given the above treatment is given by the massless bubble effect 
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and is far too small to be seen experimentally [s^ . However, it is interesting that quantum effects do predict such 
a difference, without any free parameter, modifying one of the key features of classical general relativity. Moreover, 
this phenomenon represents another demonstration that effective field techniques can make well-defined predictions 
within quantum gravity. 


9. MASSLESS PARTICLE SCATTERING 

There are also calculations of the gravitational scattering involvin g all massless particles. These include a remarkable 
calculation of graviton-graviton scattering by Dunbar and Norridge[^, a calculation of the scattering of a real scalar 
by the same authors [13 and the scattering of two non-identical scalars [blj obtained using the methods of [iOl. The 
character of the results are similar to the reactions discussed above, so we can refer the reader to the original papers 
for the results. However, we do want to comment on a couple of features of the massless amplitudes. 

One interesting feature is that there are no “classical” corrections in the massless amplitudes. We have seen this 
partially in the previous section, where the square-root non-analyticity is associated with the massive leg only not 
with the massless field. Amplitudes of totally massless particle only involve logarithms. Therefore the loop expansion 
is here strictly the h expansion, and one does not build up classical gravitational solutions in the intermediate states. 
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A second interesting feature is the ubiquity of infrared singularities and non-local effects. All of these scattering 
amplitudes can be decomposed into coefficients times the bubble, triangle and box diagrams. For massless particles 
these are particularly simple: 
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Of the many divergences that are displayed in these equations, only the 1 /e in the bubble diagram I 2 (s) is of ultraviolet 
origin. Therefore only it and the related +2 in the same amplitude represent local physics. All the other 1/e term 
in the other amplitudes represent infrared divergences, and the logarithms represent non-local effects. This tells 
us that almost all loop processes in gravity are either infrared divergent or non-local. While these features are 
well understood in scattering amplitudes, they are far from well understood in other gravitational settings such as 
cosmology and classical solutions. The IR portions of loops point to new phenomena in gravitational physics. 


10. CONCLUSION 

Although a renormalizable theory which merges general relativity and quantum mechanics has yet to be identified 
by experiment, we have shown above that when treated as a nonrenormalizable effective field theory, quantum gravity 
is a very successful theory. An EFT represents an expansion in powers of derivatives (energy-momentum) divided 
by a scale parameter, needed to make the expansion parameter dimensionless. In the case of chiral perturbation 
theory, which is a very successful picture of low energy QCD, the chiral scale parameter A^ ~ dTrE^. ~ 1 GeV so 
that xpt is valid at energies ~< 500 MeV. On the other hand, the scale parameter in the case of quantum gravity is 
hgrav ~ Mpi ~ 10^® GeV so that effective gravity is valid at any energy which is reachable in current accelerators. 
This means that higher order gravitational counterterms have essentially no influence on current experiments. On the 

/ _ 2 _ 2 

other hand, the existence of gravitational loop effects means that various non-analytic terms such as y or log 
are present in transition amplitudes, indicating, after Fourier transform, the presence of long distance—l/r" with 
n = 2, 3,4...—effects in quantities such as the energy-momentum tensor, the metric tensor, the interaction potential, 
etc. Despite the fact that these are loop effects, they are expected from classical physics arguments and in the 
case of and these forms agree with well known classical solutions. In addition loop effects produce quantum 
mechanical corrections to the classical results of order h/ (TOr)"+^ with n = 2, 3, A.... Indeed, it is the nonlocal quantum 
effects which are the purest manifestation of the EFT. 

We have evaluated a series of scattering amplitudes and focused on the parts of the calculation that the effective 
theory is capable of calculating. We found in each case that the quantum corrected amplitude is well defined when 
gravity is treated as an effective field theory. Some of the results that we displayed were: 

i) Gravitational corrections to the energy-momentum tensor of a massive system: one loop gravitational cor¬ 
rections were calculated for matrix elements of the energy-momentum tensor < P2|2/ii/(a:)|pi >. The results 
were found to agree with the classical forms of the gravitational energy-momentum tensor for both spinless 
and spin 1/2 systems. Gonverting to the metric tensor by use of the (linear) Einstein equation, the spinless 
case was found to agree with the Schwarzschild solution, while in the case of spin 1/2 the corrections were 
shown to match the Kerr solution. Quantum corrections to these classical results were determined in both cases. 

ii) Gravitational scattering of two massive particles: one-loop gravitational corrections to the potential which 
characterizes the interaction of particles with mass mi and m 2 were evaluated. At lowest order the potential 
is simply the classic Newtonian result. However, at higher order modifications are found and once again there 
exist both long distance classical and quantum mechanical corrections. Here we found a universal soft theorem, 
such that the form of the quantum correction is universal, coming as a reflection of the soft theorems of tree 
amplitudes. 
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iii) Gravitational scattering of a massive and massless particle: one loop gravitational corrections were calculated to 
the interaction of a massive test particle with a massless scalar and photon. In the case of a Feynman diagram 
calculation involving a massless system the higher order energy-momentum tensor is not well defined and the 
calculation of the scattering amplitude requires the use of both photon and graviton mass regulators. Individual 
diagrams have numerous divergences but when the total is calculated, there exists a subtle cancelation of the 
various divergences so that the final result is finite. However, modern on-shell helicity amplitude methods were 
shown to provide a significantly simplified route to this result. Here we found that massless particles no longer 
follow null geodesics, and that different types of massless particles have different trajectories. Both of these 
results are deviations from classical behavior. 

iv) Although we have not highlighted this in the discussion above, one can also see from these results that there is 
not a form of a “running” coupling G{E) in the effective theory [hlj. There is no universality of the quantum 
corrections which could have been absorbed into a running coupling. This result is totally expected in the 
effective field theory. However it is worth stating, as there are many attempts in the literature to define such a 
running G. 

Of course, general relativity is not primarily concerned with scattering amplitudes. We have started with these 
because scattering is what perturbative quantum field theory does best. However, the next challenge becomes to extend 
these effective field theory techniques to other solutions of general relativity. The metric calculations described above 
are a start down this path. There have also been applications to cosmology [h^. But much remains to be done to 
understand the low energy quantum predictions that can be calculated in the effective field theory. 
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Appendix A 

In this section we derive the various couplings to be used in our calculation. We begin with the scalar field, whose 
matter action is 


= y (137) 

where is the covariant derivative with respect to the background field. (In our case quantizing about flat space 
we have = 9^-) 

The gravitational coupling of a spin-0 particle is found by expanding the minimally coupled scalar field matter 
Lagrangian 


(138) 

in terms of the gravitational field which is a small fluctuation of the metric about flat Minkowski space defined as 

T (139) 

with K = V327rG oc 1/Mp. The inclusion of this factor k in the definition of the graviton field gives this field a 
mass-dimension of unity and thus yields a kinetic term of standard normalization without a dimensionful parameter. 
For matter interactions, this choice is convenient since the order of k keeps track of the number of gravitons involved 
in an interaction. Once the action is written in terms of the expansion of the graviton field, all indices are understood 
to be lowered or raised using the Minkowski metric We also require the expansion of the inverse metric and 
square root of the determinant of the metric tensor— 

=1 + + y + 0{n^). 


(140) 
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Then, expanding in powers of k, we find 
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where represents the trace so the one- and two-graviton vertices are identified as 
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FIG. 5: The one- and two-graviton couplings 
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FIG. 6: The three graviton vertex 
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We also require the energy-momentum tensor for gravitons which leads to the triple-graviton vertex[ 
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+ [q\q^{'napl^''''is + q-tsi^'''ap) + qxq'' + v-fsi^^^a/s) 

q ’75 ct/ 9 ) q^ q q ~k 

+ [2qYr''’^^YsMk - qT + i^^’cpYs.xaY - qf 

— ~fSlaP,Xak^ — 

+ q^{I'^^'o,pYs,Y + /a/3.a"/"'^ 75 ) + q^q.^o^PiS + Y5,XpI''"'o,p)] 

+ [{Y + {k- qf) (^r^^’o,pYs,Y + r^'o,pYs,a^ - \v^''PcPn& 

- {Yq^sl‘''''ap + {k - qfqapl^'''^s)]] (144) 

For the case of spin-1/2 we require some additional formalism in order to extract the gravitational couplings, which 

is necessary because the Dirac algebra { 7 “, 7 ^} = 2q°‘^ is defined with respect to the Minkowski flat space metric. In 
this case the Dirac matter Lagrangian coupled to gravity reads 


x/^Cm = x/^i^ 


■e^a{l°',Df,} - m 


ill 


(145) 


and involves the vierbein a which links global coordinates with those in a locally flat space. The vierbein is in some 
sense the “square root” of the metric tensor and satisfies the relations 


qah — qfJ:!^ 
e^eYg^’' = qab 


e^aeYq^^ = 9^'' 
e>^aeYgp. = q^'^. 


The covariant derivative is 


— 2 W ^ ^ & "“C ^ 


with cr‘=^ = 


c d 

7,7 


and the partial derivative acts only on spinors and in such a way that 


= tjj df_,ip - {df,Y) '0- 


Putting everything together, we find then 
x/^Cm = x/^i’ 


1 - I eY' u^Yb Vac { 7 “', ^ 


' 0 - 


(146) 

(147) 

(148) 

(149) 


The spin connection ojYi, Vac can be derived in terms of vierbeins by requiring 0^6^°“ = 0 and by antisymmetrization 
in /r O iz in order to get rid of Christoffel symbols^. The result is: 

ujYbqac={^eY{d^eY-d,eY) + '^eYeYe/dpeY) - (fc O c) (150) 

In order to derive the Feynman rules we expand the ingredients in Eq. (I149p that contain graviton couplings, that is 
we need e^a and 10 ^ 1 , Vac expanded up to 0{Y) 

= K-1 ^ + ... 

uJ^bVac = '^dbh^!^} + y - '\hP'>Ydph^^} + ^- (b ^ 

(151) 


^ For our purposes we shall use only the symmetric component of the vierbein matrices, since these are physical and can be connected 
to the metric tensor, while their antisymmetric components are associated with freedom of homogeneous transformations of the local 
Lorentz frames and do not contribute to nonanalvticitv fbSl 
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After these expansions are employed, we no longer need to distinguish between Latin Lorentz ind ices and Greek 
covariant indices and can use the Minkowski metric to lower and raise all indices. 

The matter Lagrangian then has the expansion—(note here that our conventions are 75 = and £^123 _ 

+ 1 ) 


and the corresponding one- and two-graviton vertices are found to be 


^ iP2: Pit 


-M 

mh^ 

m)=iK^ 


+P2)v+lv{pi+P2)^-'npu\^i^l+l/)2)- 


1 

16 


- \ Q(|^1 + 1^2) - rri^ Ppu,pa 

Vpi'{lp(pi + P2)a + 'Jcripi +P 2 ) 

+ Vpa{lp{Pl +P2)u +lu(j>l +P2)p^ 


+ Y^(pi +P2ypHk^ ,fiip I* 

e,p (7 “t“ 

~f TAT5 ^pvpa T 


(152) 


(153) 


Appendix B 

Above we have argued that the scattering amplitude which is ultimately related to observables in quantum field theory 
is a physical quantity while the potential we have given is not and depends on the gauge, the choice of coordinates, 
and on the way the iteration is performed, be., on the way we do the matching. While the classical component of our 
potential is in fact plagued by these ambiguities, the quantum part is unique since it is unaffected by how we perform 
the matching and a quantum field theory calculation in any gauge would result in the same result |3l|. 

In this appendix we demonstrate how we can recover the classical equations of motion from our scattering amplitudes 
by setting up the Einstein-Infeld-Hoffmann (EIH) Lagrangianf^. The EIH Lagrangian is itself dependent on the 
choice of coordinates, but can be expressed in the center of mass frame {P = pi = —ps, r = ri — r^) in a general 
way as[^ 


Leih = T —V 


where the kinetic energy to NLO in the nonrelativistic expansion reads 

p2 p2 pA p4: 

T = _^_ 

2mi 27712 8mf Sm^ 


(154) 


(155) 


and the potential is 


V = V^^'> + 


( 156 ) 
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with 


1/(1) 


1/(2) 


GtoiTO2 

r 



1 

2 



TOi77l2 miTO2 


(mi + m2)2 


77111712 


jP ■ ff 

17111712 


(1 - 2a) 


G^TOim2(mi + m 2 ) 

^ ■ 


(157) 

(158) 


The parameter a parameterizes the choice of coordinates used, with a = 0 being the gauge of the original EIH result. 
The coordinate change 


r —7> r 


^ ^ G(mi + m 2 ) ^ 


(159) 


which implies 


P —7 P + a 


G(mi + m 2 ) 


[P-{P- f) f] 


(160) 


then brings the original EIH Lagrangian into the form above, which is the most general result. 

Since we perform our matching on-shell, i.e., we use the on-shell one-graviton exchange amplitude to define the 
leading order 0{G) potential, terms proportional to P ■ f never arise, meaning that our result corresponds to a gauge 
such that the coefficient of the structure 


jP-f)" 

mim2 


in Eq. (115711 vanishes. That is the case if and only if the gauge parameter is 


a = — 


mim2 


2(mi -I- m2)2 


whereby the EIH potential becomes 

1/(1) = 
H(2) = 


Gmim2 


1-f 


^ mim2 
(mi -I- m 2 )^. 


3 (mi -I- m 2 )^ 
mim2 
G^mim2(mi -I- m 2 ) 
^ ■ 


^+2 


(161) 


(162) 

(163) 


Comparing the EIH potential in this gauge of Eq. (11621) with the long distance component of the leading order 
spin-independent potential in Eq. m we find full agreement for the relativistic corrections to the 0{G) potential. 
However, comparing the EIH potential 17(^) in this gauge—Eq. (11631) — with the classical component of our spin- 
independent potential in Eq. dMl) we see that the two do not agree! The reason for this discrepancy is that we elected 
to use nonrelativistic forms when we performed the second Born iteration of the leading order potential in Eq. (1981) . 
This procedure, however, is not self-consistent when we are interested in equations of motion at NLO, and we must 
account for the leading relativistic corrections in performing the iteration. In particular, we must use expressions for 
the potential and the propagator in Eq. (1981) which include these leading relativistic terms^ 


Pf 


Oy(l) 

^NLO 


Pi ) 


47rGmim2 


1 + 


Pi 


2mim2 


Pf 3(mi -bm2)2 


2mim2 


(164) 




Pa 

2mr 


J2- 

2mr 


le 


i-f 


f^ + JL 

\ 4m2 4m2 


1 - 3- 


mim2 


( 165 ) 


Tlie subscript NLO in this sections refers to the iteration being performed at NLO in the relativistic expansion. 
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which yields a second Born iteration amplitude 

d^£ dTrGm 17712 




47rGmim2 


(27r)3 |p/-£|2 

2mr 2mr 




1 + 


(PO + £'^) ( I , 7 77717772 ^ 


777 i7772 \4 ' 4 777^ / 

L 777r G^rri^m 


■ A /^2 2 2 "''I 

247rG 77111712 --^ 


Po mi + m2 


1"C2 I ^ 7(7771 + 7772)^ j ^ 


mim2 


(166) 


Subtracting this iterated amplitude, which includes all corrections to NLO, from the scattering amplitude (9) 
of Eq. (iroi) we find the second order potential 


(r) — — f ^ ^ 

nloT')- J (2^)3e 


0 t /( 2 ) 

y T 


- “Amp)^)^o(g) 

mim2 


.(2) 


g-zq.r- g2^17772 


(27r) 


(7771 + 7772 ) + 


7771 + W2 1 5 


= 1 + 


77717772 \ G^777 i 7772 (tTTi + 7772) 41G^777 i7772/7 


(ttt-i + 7772 )^ 


2r2 


IOttc^ 


(167) 


and observe that now the classical component agrees with the G(G^) EIH potential of Eq. (I163E 

Thus we have shown that if we consistently take into account the and Gmjr corrections beyond Newtonian 
physics we reproduce the EIH Lagrangian in a certain gauge. Erom the resulting EIH Lagrangian we could evaluate 
observables such as the precession of the perihelion of Mercury which must clearly be independent of gauge. The 
inclusion of the corrections is required since the equations of motion can be used to describe bound states where 
~ Gmjr by the virial theorem. 

However, our methods are clearly clumsy for the calculation of classical observables. Recently, Goldberger and 
Rothstein have developed an effective field th eory of gravity which is optimized for calculating classical observables 
of bound states called NRGR[6^,[63,[il|,[i^,[Zfl) Zl| Here the external particles are static sources so that no loops 
are to be calculated in their theory when calculating classical observables since the only propagating particles present 
are gravitons which are massless and thus the loop expansion in NRGR corresponds to an expansion in h. In the 
NRGR framework the ^in-dependent classical equations of motion were calculated recently to NLO by Porto and 
Rothstein so that we will not continue here to evaluate the corresponding spin-dependent classical 

potentials consistently taking into account all relativistic 0{v^) effects in the iteration. 
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